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Specific  rate  constants,  k(E,  J),  and  thermal  capture  rate  constants,  kaJT),  are  determined 
by  statistical  adiabatic  channel  model/classical  trajectory  (SACM/CT)  calculations  for 
unimolecular  dissociation  and  the  reverse  association  reactions  of  representative  polyatomic 
molecular  ions.  Simple  short-range  valence/long-range  ion-induced  dipole  model  potentials 
without  reverse  barriers  have  been  employed,  using  the  reactions  C8H10+  C7II7+  +  CH7 
and  C9H,2+<s>C7H7+  +  C2H5  as  illustrative  examples.  Simplified  representations  of  k(E) 
and  kcap(T)  from  rigid  activated  complex  Rice-Ramsperger-Kassel-Marcus  (RRKM) 
theory  are  compared  with  the  SACM/CT  treatment  and  with  experimental  results.  The 
Massey  parameters  of  the  transitional  mode  dynamics,  for  the  systems  considered,  are 
smaller  than  unity  such  that  their  dynamics  is  nonadiabatic  while  the  dynamics  of  the 
conserved  modes  is  adiabatic.  Because  of  the  long-range/short-range  switching  character 
of  the  potential,  simple  rigid  activated  complex  RRKM  theory  cannot  be  used  without 
modifications.  The  effects  of  a  shifting  of  the  effective  bottle-neck  of  the  dynamics  with 
increasing  energy  towards  smaller  interfragment  distances  in  the  present  cases  are  ampli¬ 
fied  by  a  shift  into  a  range  of  increasing  anisotropy  of  the  potential.  As  a  consequence, 
the  thermal  capture  rate  constants  markedly  decrease  with  increasing  temperature. 
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1.  Introduction 

In  the  previous  articles  of  this  series  [1-6],  we  have  treated  unimolecular  bond 
fission  and  the  reverse  association  reactions  by  a  combination  of  the  statisti¬ 
cal  adiabatic  channel  model  (SACM)  and  classical  trajectory  (CT)  calculations. 
The  dissociative  separation  of  two  species  and  the  reverse  associative  capture 
process  forming  a  bound  adduct  were  followed,  in  reduced  dimensionality,  by 
explicitly  treating  the  dynamics  of  the  transitional  modes  (external  rotations  of 
the  fragments)  while  conserved  modes  (internal  modes  of  the  fragments)  were 
handled  adiabatically.  The  SACM/CT  approach  is  able  to  establish  a  direct  link 
between  the  potential  energy  surfaces  (PESs)  and  thermal  rate  constants  for 
capture  kap(T)  (or  the  corresponding  high  pressure  rate  constants  for  recombi¬ 
nation  and  the  reverse  dissociation)  such  as  demonstrated  in  [1-5].  In  addition, 
specific  rate  constants  k(E,  J )  for  unimolecular  bond  fissions  as  a  function  of 
the  energy  E  and  the  angular  momentum  (quantum  number  J)  can  be  derived 
in  a  way  which  is  consistent  with  the  thermal  capture  rate  constants,  see  [6]. 

Calculating  the  ./-dependence  of  k(E,  J )  is  an  important  element  of  the 
analysis,  since  thermal  capture  involves  broad  distributions  over  E  and  J.  In 
addition,  complex-forming  bimolecular  reactions,  to  which  the  same  approach 
maybe  applied,  involve  subtle  competitions  between  E-  and  J-dependent  com¬ 
peting  association  and  dissociation  processes  of  the  adducts  [7].  Generally,  the 
./-dependence  is  handled  in  a  more  simplified  manner  than  presented  here  and 
the  employed  average  /-treatments  are  not  always  adequate,  see  e.g.  [8]. 

It  is  obvious  that  the  kinetic  properties  are  governed  by  PESs  such  that 
direct  links  between  potentials  and  rate  parameters  need  to  be  established  in 
a  transparent  way.  For  this  reason,  we  have  calculated  thermal  capture  rate  con¬ 
stants  [1-5]  and  the  corresponding  specific  rate  constants  [6]  for  a  series  of 
interfragment  potentials  such  as  ion-dipole  [1],  ion-quadrupole  [2],  dipole- 
dipole  [3],  atom-linear  rotor  model  valence  [4],  and  linear  rotor-linear  rotor 
model  valence  [5]  potentials.  By  using  such  model  potentials,  insight  into  what 
specifically  controls  the  reactivity  is  developed.  SACM-  and  CT-treatments  of 
the  transitional  mode  dynamics  were  shown  to  coincide  in  the  limit  of  adia¬ 
batic  and  classical  behavior.  Nonadiabatic  dynamics  and  quantum  effects  were 
also  investigated.  It  was  shown  [9]  that  the  SACM/CT  concept  works  down  to 
temperatures  in  the  milli-  and  even  micro-Kelvin  range,  particularly  if  Coriolis 
coupling  effects  are  properly  accounted  for  [10].  The  presence  of  open  elec¬ 
tronic  shells  also  becomes  relevant  at  low  temperatures  and  needs  to  be  taken 
into  account,  see  e.g.  [11]  and  [12], 

The  transitional  mode  potentials  considered  in  [1-6]  were  partly  of 
“mixed”  (ion-permanent  plus  induced  dipole  [1],  ion-quadrupole  plus  in¬ 
duced  dipole  [2],  dipole-dipole  plus  dispersion  [3]),  partly  of  “homogeneous” 
valence  character  [4, 5].  The  former  potentials  [1-3]  were  characterized  by  re¬ 
gions  of  differing  anisotropies  while  the  latter  [4, 5]  had  uniform  anisotropies. 
This  raises  the  question  which  parts  and  properties  of  the  potential  are  most 
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relevant  for  the  dynamics.  Answers  are  provided  by  inspecting  the  positions 
of  maxima  of  adiabatic  channel  potential  curves  in  SACM  [13]  or  positions 
of  transition  states  in  flexible  variational  transition  state  theory  (VTST)  [14]. 
With  both  approaches  one  observes  a  smooth  “transition  state  shifting”  to¬ 
wards  smaller  interfragment  distances  with  increasing  energy  such  that  regions 
of  different  anisotropy  characters  may  be  sampled.  There  may  also  be  more 
abrupt  “transition  state  switching”  if  the  PESs  have  small  barriers  or  reefs 
along  the  reaction  coordinate  such  as  found  in  ab  initio  calculations  of  PESs 
for  small  molecules  like  H02  [15],  HOCO  [16],  or  H202  [17].  Whether  tran¬ 
sition  state  switching  also  occurs  without  barriers  or  reefs,  remains  an  open 
question  [8]  and  one  would  like  to  find  an  answer  by  SACM/CT  calcula¬ 
tions  on  the  given  PESs.  However,  there  is  a  problem:  if  the  transitional  mode 
dynamics  is  nonadiabatic,  adiabatic  channel  potential  curves  loose  their  signif¬ 
icance  and  sampling  of  phase  space  volumes  in  simple  VTST  likewise  looses 
its  rigorous  relation  to  dynamics.  Even  the  most  elaborate  variable  reaction  co¬ 
ordinate  VTST  [18],  requires  “dynamical  corrections”  to  be  determined  by  CT 
calculations  [19]. 

In  the  present  work  we  consider  PESs  of  mixed  character,  i.e.  valence-ion- 
induced  dipole  potentials,  which  link  regions  of  quite  different  anisotropies. 
For  such  potentials,  a  switching  between  inner  tight  and  outer  orbiting  transi¬ 
tion  states  has  been  invoked  for  molecular  ions  [20]  when  transitional  mode 
frequencies  along  the  reaction  coordinate  were  assumed  to  change  in  a  way 
which  strongly  differs  from  the  changes  of  the  minimum  energy  path  (MEP) 
potential.  The  question  remains  [8,20,21]  to  what  extent  true  transition  state 
switching  is  a  consequence  of  the  chosen  potential  and  of  the  general  approach, 
and  whether  transition  state  switching  is  real  for  molecular  ions  besides  the 
normal  transition  state  shifting.  Because  of  the  general  importance  of  this  ques¬ 
tion  we  apply  the  SACM/CT  approach  to  polyatomic  ions  with  model  PESs  of 
the  described  mixed  character.  We  focus  attention  on  an  internally  consistent 
calculation  of  specific  rate  constants  k(E,  J)  and  thermal  capture  rate  con¬ 
stants,  kcav(T).  In  a  separate  publication  [22],  we  also  provide  an  SACM/CT 
calculation  of  product  energy  distributions  (PEDs)  which  are  internally  consis¬ 
tent  with  the  present  calculations.  Our  results  for  k{E,  J),  kcap(T)  and  PEDs 
show  all  signs  of  a  transition  state  shifting  or  transition  state  switching  but 
cannot  distinguish  between  them.  With  increasing  E,  J,  or  T,  dynamical  con¬ 
straints  caused  by  the  increasing  anisotropy  of  the  PES  at  shorter  interfragment 
distances  lead  to  results  which  would  correspond  to  increasingly  rigid  transi¬ 
tion  states.  Since  the  chosen  reaction  systems  show  nonadiabatic  transitional 
mode  dynamics,  however,  for  principal  reasons  it  would  be  difficult  if  not  im¬ 
possible  to  rigorously  define  transition  states.  Since  there  is  also  no  real  need  to 
look  for  them,  the  question  of  transition  state  shifting  or  switching  finally  be¬ 
comes  irrelevant.  Previous  transition  state  switching  treatments  [20]  have  not 
provided  such  a  consistent  approach  as  the  full  ./-dependence  was  not  elabo¬ 
rated. 
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The  outline  of  the  present  article  is  as  follows.  First,  we  construct  reduced 
dimensionality  model  PESs  for  transitional  modes  in  molecular  ion  dissoci¬ 
ations.  We  use  the  dissociations  of  alkylbenzene  ions  as  practical  examples 
because  conventional  modelling  of  k(E )  for  these  systems  by  rigid  activated 
complex  Rice-Ramsperger-Kassel-Marcus  (RRKM)  theory  and  experimental 
results  are  available.  Classical  trajectory  calculations  are  performed  on  these 
PESs,  leading  to  k(E,  J )  and  kc.^,(T)  (and  PEDs(E,  J),  see  [22]).  A  minimum 
of  fitting  of  the  potentials  is  necessary  to  bring  modelled  k(E,  J )  into  agree¬ 
ment  with  the  experiments.  The  corresponding  rate  constants  kC!tp(T )  then  are 
obtained  by  thermal  averaging  of  k(E,  J )  and  again  compared  with  experi¬ 
ments.  Various  fitting  procedures  of  both  the  potential  and  the  rate  constants 
are  inspected  and  compared  with  traditional  RRKM  modelling  with  the  goal 
to  identify  the  most  economic,  but  still  realistic  treatment  in  practical  applica¬ 
tions,  in  particular  with  respect  to  the  relation  between  specific  and  thermal  rate 
constants. 


2.  Model  potentials  for  the  dissociation  of  polyatomic 
molecular  ions 

The  SACM/CT  approach  treats  the  dynamics  of  the  transitional  modes  explic¬ 
itly  such  that  an  analytical,  reduced  dimensionality,  PES  in  these  coordinates 
is  required.  The  potential  of  the  conserved  modes  is  simplified  assuming  that 
these  mostly  vibrational  modes  are  adiabatic  and  that  their  frequencies  are 
close  to  the  fragment  frequencies.  The  latter  approximation  could  be  relaxed 
by  correcting  for  the  change  of  the  frequencies  along  the  interfragment  distance 
and  choosing  frequencies,  e.g.,  near  to  the  positions  of  the  centrifugal  barriers. 

The  SACM/CT  method  calculates  trajectories  in  the  transitional  mode  po¬ 
tential.  Initial  conditions  for  these  CT  calculations  are  two  fragments  at  large 
separation.  Combining  their  individual  rotational  angular  momenta  into  one 
total  angular  momentum,  the  number  W0(E,  J )  of  their  quantum  states  of  total 
energy  E  and  total  angular  momentum  (quantum  number  J )  is  calculated  first. 
During  the  approach  of  the  two  fragments,  “capture”  occurs  when  the  inner 
turning  point  of  the  relative  motion  is  clearly  inside  the  well  of  the  bound 
adduct. 

Neglecting  the  anisotropy  of  the  potential,  a  capture  probability  wrST(E,  J ) 
then  is  defined  which  denotes  that  fraction  of  the  W0(E,  J)  quantum  states 
which  lead  to  capture  in  the  artificially  isotropic  potential.  In  general,  the 
product  JVFST(E,  J )  =  vjvst(E,  J )  W0(E,  J )  will  be  smaller  than  W0(E,  J)  be¬ 
cause  dynamical  constraints  arise  in  the  form  of  centrifugal  barriers,  E0(J). 
W?ST(E,  J )  is  termed  the  “number  of  open  channels”  (or  of  “activated  complex 
states”)  in  phase  space  theory  (PST).  One  should  note  that  our  version  of  PST 
differs  from  other  versions  of  PST  [8, 23],  e.g.  from  the  simplest  form  in  which 
centrifugal  barriers  are  completely  neglected.  Our  version  of  PST  also  differs 
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from  orbiting  transition  state  PST  (OTS/PST)  where  centrifugal  barriers  E0(J) 
are  calculated  with  the  long-range  electrostatic  potential  only  [24],  The  present 
PST  provides  a  rigorous  upper  limit  of  the  capture  rate  which  would  not  be 
the  case  in  OTS/PST  which  neglects  short-range  components  of  the  PES,  see 
below. 

Including  the  anisotropy  of  the  real  potential  introduces  additional  dy¬ 
namical  constraints  such  that  the  product  of  the  corresponding  capture  prob¬ 
ability  w(E,  J)  and  W0(E,  J),  i.e.  W(E,  J)  =  w(E,  J)W0(E,  J),  is  smaller 
than  fVPST(E,  J).  Convoluting  the  numbers  of  open  channels  W(E,  J)  of  the 
transitional  modes  with  the  contributions  of  conserved  modes  completes  the 
calculation  of  the  total  number  of  open  channels  Wtot(E,  J).  Thermal  aver¬ 
ages  over  Wtot(E,  J)  determine  the  thermal  capture  rate  constants,  kcap(T).  At 
the  same  time,  Wtot (E,  J )  also  enters  the  expression  for  specific  rate  constants 
k(E,  J)  =  JVm(E,  J)/hp(E ,  J)  such  as  employed  in  statistical  unimolecular 
rate  theory  [8],  p(E,  J)  in  the  usual  way  denotes  the  rovibrational  density  of 
states. 

SACM/CT  calculations  in  the  past  have  been  performed  on  ab  initio  po¬ 
tentials,  see  e.g.  [15]  and  [17],  and  on  model  potentials,  see  [4]  and  [5],  In  the 
present  work,  we  are  interested  in  polyatomic  ions  for  which  ab  initio  poten¬ 
tials  are  available  only  in  rare  cases  (see,  e.g.,  [25-27]).  For  this  reason,  we 
design  simple  model  potentials  which  connect  short-range  valence  with  long- 
range  electrostatic  potentials.  The  transition  between  the  two  components  of 
the  PES  can  be  smooth  like  calculated  for  LiH2+  dissociation  [26]  or  abrupt 
like  derived  for  C6H6+  dissociation  [25],  Unfortunately,  often  ab  initio  calcula¬ 
tions  of  PESs  for  molecular  ions  (like  the  PES  for  N4+  dissociation  from  [27]) 
are  not  extended  to  that  range  of  interfragment  distances  where  the  transition 
from  short-range  to  long-range  potential  occurs  such  that  there  is  no  a  priori 
guide  to  the  design  of  the  model  potential.  For  the  C8H10+  and  C9Hi2+-systems 
of  primary  interest  in  the  present  article,  we  decided  to  construct  a  model  po¬ 
tential  with  smooth  short-range/ long-range  switching.  Because  of  the  model 
character  of  the  potential,  we  had  to  allow  for  the  introduction  of  one  fitting 
parameter  which  was  chosen  to  be  an  anisotropy  amplitude  factor.  However, 
this  fitting  parameter  was  fixed  by  means  of  a  single  measurement  and  no  more 
fitting  was  allowed  for  comparison  with  other  measurements. 


2.1  Radial  part  of  the  potential 

It  is  important  to  realize  which  property  of  the  rate  constants  is  governed  by  the 
radial  part  of  the  potential  and  which  is  connected  to  the  anisotropy.  For  this 
reason,  we  first  consider  the  radial  part  of  our  model  potential  only.  The  short- 
range  component  of  this  potential  is  assumed  to  be  of  Morse  type 


ySR(r)  =  D{  exp[— 2/3(r  —  re)]  —  2  exp[— /3(r  —  re)]} . 


(2.1) 
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In  the  present  work  we  consider  a  long-range  induced-dipole  component  of 
the  potential  of  the  form 

VLR(r )  =  - aq2/2r 4  (2.2) 

(for  our  nomenclature  see  [1-6]),  but  other  long-range  potentials  could  be 
treated  in  an  analogous  way.  Including  the  ion-dipole  term  has  little  effect  on 
the  present  calculations  since  the  polarizability  is  large  and  the  dipole  moment 
is  small.  We  assume  that  the  center-of-mass  (c.o.m.)  distance  r  between  the 
fragments  can  serve  as  the  reaction  coordinate  and  that  the  radial  potential  F(r) 
is  the  MEP  potential. 

There  are  various  ways  to  represent  the  switching  of  V(r)  between  FSR(r) 
and  FLR(r).  For  example,  one  may  try  to  use  an  expression  [26] 

V{r )  «  VSR(r)F(r)  +  FLR(r)[  1  -  F(r)]  (2.3) 


where  F(r )  is  an  exponentially  decaying  switching  function.  In  the  present 
work  a  more  sophisticated  switching  procedure  was  employed  for  reasons 
which  become  clear  below.  We  modified  VLR(r)  in  two  stages.  First,  Eq.  (2.2) 
was  replaced  by 


VLR(r)  =  -aq2/2r* 


(2.4) 


where  rB,  instead  of  the  c.o.m.  distance  r,  denotes  the  length  of  the  breaking  (or 
of  the  forming)  bond.  We  then  removed  the  discontinuity  of  the  denominator  of 
Eq.  (2.4)  by  adding  a  term  ocq2/2D,  such  that 


FLR(r)  «  - 


aq2/2 

(r  —  ArB)*  +aq2/2D 


(2.5) 


where  ArB  =re  —  rBe  (the  subscript  e  denotes  equilibrium  values  of  r  and  rB). 
Second,  instead  of  an  exponential  switching  function,  we  used 


V(r)  =  FSR(r)[l  -  kF]  +  FLR(r)F 


(2.6) 


with  FSR(r)  from  Eq.  (2.1),  and  VLR(r)  from  Eq.  (2.5).  F  and  k  are  defined  as 


F=  l/{l+exp[-2^(r-rsw)]}, 


(2.7) 


and 

^  =  1 - ^ - 4  .  (2.8) 

4£>(rsw-ArB)4 

rsw  denotes  that  value  of  r  where  FLR(r)  from  Eq.  (2.5)  intersects  with 
FSR(r)  from  Eq.  (2.1).  The  radial  MEP  potential  from  Eqs.  (2.6)-(2.8)  was  de¬ 
signed  to  produce  a  smooth  monotonic  transition  from  Eq.  (2.1)  to  Eq.  (2.2). 
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Fig.  1.  Radial  potential  V(r)  for  the  reaction  C8H10+  C7H7+  +  CH3  (short-range  Morse 

potential  FSR(r)  from  Eq.  (2.1);  modified  long-range  ion-induced  dipole  potential  VLR(r) 
from  Eq.  (2.5);  model  potential  V(r)  from  Eqs.  (2.6)— (2.8);  maximum  angular  poten¬ 
tial  f'anisoM  =  Cx  exp[-^(r  -re)]  with  CJD  14  from  Eq.  (2.15)  for  ion-dipole  type 
anisotropy;  for  molecular  parameters,  see  Appendix,  all  quantities  represented  as  V(r)/D). 


In  order  to  achieve  the  smooth  transition  for  various  systems,  the  expres¬ 
sion  for  k  needed  some  flexibility;  e.g.  Eq.  (2.8)  was  found  to  be  adequate 
for  the  reaction  C8Hi0+  — >  C7H7+  +  CH3,  whereas  Eq.  (2.9)  was  preferred  for 
C9H12+  -»  C7H7+  +  C2H5, 

^  =  1 - (2.9) 

D  frsw  —  Ar B) 

because  it  gave  a  smoother  transition  than  Eq.  (2.8).  We  illustrate  the  modelled 
MEP  potentials  for  the  dissociations  of  C8H10+  and  C9H,2+  in  Figs.  1  and  2. 
We  use  the  molecular  parameters  specified  in  Appendix  A.  The  modelled  MEP 
potential  shows  smooth,  but  relatively  fast,  transitions  between  VSR(r)  from 
Eq.  (2.1)  and  VLR(r)  from  Eq.  (2.5).  It  gives  major  deviations  from  the  short- 
range  valence  potential  above  about  V{r)/D  ~  —0.5,  similar  to  the  ab  initio 
PES  for  the  LiH2+  -»•  Li+  +  H2  reaction  [26]  which  shows  major  deviations 
from  a  Morse-type  valence  potential  at  V{r)/D  above  about  —0.5.  However,  in 
the  case  of  the  production  of  an  atomic  ion,  the  long-range  potential  at  V(r)  / D 
~  —0.5  does  not  intersect,  but  only  comes  close  to  the  short-range  potential. 
For  the  considered  polyatomic  ions,  long-range  potentials  from  Eq.  (2.2)  would 
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Fig.  2.  As  Fig.  1,  for  C9H12+  <3>  C7H7+  +  C2H5  (with  C,/Z)  «  12). 


intersect  the  short-range  potential  at  V{r)/D  above  about  —0.05.  For  this  rea¬ 
son,  replacing  r  by  rB  in  Eq.  (2.2)  appears  more  physical,  leading  to  the  results 
shown  in  Figs.  1  and  2. 

The  radial  potentials  are  all  that  is  required  for  a  PST  calculation  of 
W(E,  J).  In  particular,  one  may  directly  determine  the  centrifugal  barriers 
E0(J ),  i.e.  the  maxima  of  the  centrifugal  potentials  V(r)  +  L (L  +  l)h2 /2/xr2 
(L  —  quantum  number  of  orbital  angular  momentum).  In  the  present  case,  one 
may  apply  the  simplification  J  ~  L,  see  below.  For  the  C8H10+  system,  one 
obtains 

E0(J)/hc  «  3.18  x  10~6  cm"1  [j4/(l  +0.0169/  +  2.65  x  10"4  J179)] . 

(2.10) 

The  corresponding  expression  for  the  C9H]2+  system  is  E0(  J)/hc  ~  5.90  x 
10-7cm-'[y4/(l+9.73x  10-3y,  075  +  5.21  x  10~6./2-166)].  For  small  J,  in  both 
cases  E0(J )  of  the  ion-induced  dipole  potential  is  approached,  such  as  given  by 

E0(J)  =  [J(J+  \)h2 /2[E]2 /2aq2 .  (2.1 1) 

For  large  J,  E0(J )  takes  the  form  typical  for  Morse  valence  potentials  [28], 
i.e. 


Eo(J)  ~  CV[J(J  +  1)] 


(2.12) 
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with  v  being  in  the  range  1  to  1.5.  Eq.  (2.10)  shows  a  smooth  transition  be¬ 
tween  Eqs.  (2.11)  and  (2.12)  for  the  C8H]0+  and  C9H12+  systems.  For  the 
model  potentials  of  the  present  work,  there  is  no  sign  of  centrifugal  barrier 
switching  such  as  found  for  potentials  with  small  barriers  or  reefs,  see  e.g.  the 
H202  -*■  20H  reaction  [17]. 

Centrifugal  pseudo-partition  functions 


2ce„,  =  J2(2J+  *>  exP  [~Eo(J)/kT]  (2.13) 

j= o 


with  Eq.  (2.10)  are  easily  calculated  and  capture  rate  constants  in  PST  are  ac¬ 
cessible  through 


(2.14) 


A  special  case  of  Eq.  (2.14)  would  be  ion-induced  dipole  capture,  for 
which  Eqs.  (2. 1 1)— (2. 13),  and  (2.14)  give  the  Langevin  rate  constant  kL  = 
2jrq(a/ii)'/2.  OTS/PST  and  the  present  version  of  PST  in  this  special  case 
coincide  and  lead  to  the  same  result  of  k^I  =  kh. 


2.2  Angular  part  of  the  potential 

While  the  here  considered  long-range  ion-induced  dipole  potential  is  isotropic, 
the  short-range  valence  potential  is  always  strongly  anisotropic.  As  a  conse¬ 
quence  of  this  anisotropy,  the  free  rotations  of  the  separated  fragments  change 
into  deformation  vibrations  of  the  adduct.  The  rapidity  of  this  change  is  re¬ 
lated  to  the  particular  character  of  the  PES,  i.e.  either  a  smooth  or  an  abrupt 
switching  from  the  short-range  to  the  long-range  component  of  the  potential. 
As  discussed  above,  we  here  decided  to  assume  a  smooth  switching.  How¬ 
ever,  we  are  aware  of  the  fact  that  also  more  abrupt  switchings  like  in  C6H6+ 
dissociation  [25]  are  possible. 

We  modelled  the  anisotropy  of  the  PES  by  angular  terms  like 

V(r,  angles)  «  V(r)  +  C ,  exp[-^(r-re)]{(2  +  dln  +  d2n)/4j  (2.15) 

for  an  anisotropy  of  ion-dipole  type  (d,  and  d2  are  unit  vectors  along  the  axis 
of  the  fragment  rotors  and  n  is  the  unit  vector  in  the  direction  of  the  breaking 
bond;  Cj  characterizes  the  strength  of  the  anisotropy;  one  should  note  that  in 
a  true  ion-dipole  system  n  would  be  in  the  direction  of  the  line  connecting  the 
c.o.m.s).  Alternatively,  we  considered 

V(r,  angles)  ~  V(r)  +  C2  exp[~P(r  - re)]{[2  +  d^d2  -  3(dln)(d1n)]/4} 

(2.16) 
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for  an  anisotropy  of  dipole-dipole  type  (dx  and  d2  are  unit  vectors  like  in 
Eq.  (2.15),  n  is  the  unit  vector  in  the  direction  of  the  line  connecting  the 
c.o.m.s;  C2  characterizes  the  strength  of  the  anisotropy).  Figs.  1  and  2  include 
amplitudes  of  the  anisotropy  such  as  described  by  initially  arbitrarily  chosen 
values  of  C\.  Quite  similar  to  the  example  [26]  of  LiH2+,  our  model  PESs 
show  anisotropies  which  reach  far  beyond  the  range  of  the  short-range  valence 
potential  into  the  range  of  the  long-range  ion-induced  dipole  potential. 

Performing  CT  calculations  on  the  model  potentials,  we  realized  that  the 
angular  form  (ion-dipole  or  dipole-dipole  character)  is  only  of  secondary  im¬ 
portance.  What  matters  is  the  anisotropy  amplitude  represented  by  the  values 
of  the  amplitude  factors  C,  or  C2.  This  quantity  needs  to  be  fitted  such  as  dis¬ 
cussed  below.  In  addition,  Eqs.  (2.15)  and  (2.16)  only  model  anisotropies  of 
linear  rotor  fragments  leading  to  a  linear  rotor  adduct.  In  reality,  there  will  be 
additional  anisotropy  corresponding  to  the  top  character  of  the  fragments.  In 
the  absence  of  more  information  on  this  property  of  the  PES,  we  limited  our¬ 
selves  to  studying  the  dynamics  with  potentials  given  by  Eqs.  (2.15)  or  (2.16). 


3.  Capture  probabilities  and  numbers  of  open  channels 

Before  performing  CT  calculations  of  the  transitional  mode  dynamics  on  the 
model  potentials  described  in  Sect.  2,  it  appeared  useful  to  analyze  the  charac¬ 
ter  of  the  dynamics  with  respect  to  the  effective  mass  M  and  the  corresponding 
Massey  parameter  §  =  (2 M)x /2,  see  [6].  The  Massey  parameter  corresponds  to 
the  ratio  between  the  collision  time  in  the  radial  potential  V(r)  and  the  period  of 
motion  of  the  considered  mode.  If  f  »  1,  the  motion  of  the  mode  stays  adiabat- 
ically  in  its  quantum  states  which  can  be  calculated  by  solving  the  eigenvalue 
problem  for  fixed  interfragment  distances.  This  leads  to  the  adiabatic  channel 
potential  curves  of  S ACM  [13],  If  ^  1,  these  curves  loose  their  significance, 

there  are  frequent  transitions  between  the  adiabatic  channel  potential  curves, 
and  the  dynamics  for  £  0  approaches  the  “sudden”  limit.  Conserved  modes 

generally  are  on  the  adiabatic  side  while  transitional  modes  may  be  more  or 
less  nonadiabatic.  This  is  the  reason  why  SACM/CT  concepts  [1-6]  have  been 
developed  by  us  beyond  the  original  SACM  model  [13]. 

In  the  case  of  an  isotropic  potential,  adiabatic  (£  1)  and  sudden  (£  <«C  1) 

dynamics  of  transitional  modes  lead  to  identical  capture  rates.  It  is  the 
anisotropy  of  the  potential  which  generates  marked  differences  between  adia¬ 
batic  and  sudden  dynamics  such  as  illustrated  in  [1]  for  ion-dipole  capture 
and  in  [4]  and  [5]  for  valence  potentials.  The  Massey  parameters  for  a  valence 
potential  are  independent  of  the  energy  and  are  given  by 

?*(2  nB)'l2/hf5  (3.1) 

with  the  reduced  mass  /x,  the  Morse  parameter  f3,  and  B  =  (5,  +B2)/ 2  where 
B  ,  and  B2  are  the  relevant  rotational  constants  of  the  two  fragments  (in  energy 
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units).  Based  on  the  molecular  parameters  of  Appendix  A,  we  find  the  follow¬ 
ing  quantities  for  C8Hi0+:  B\/hc  ~  0.078  cm-1  and  B2/hc  ~  9.3  cm-1  which 
leads  to  £  =  (2M)1/2  «  0.51.  Separating  the  contributions  of  the  rotors  would 
give  £,  ~  0.066  and  §2  ~  0.72.  For  C9H12+  with  B2/hc  «  0.73  cm-1,  we  ob¬ 
tain  £j  ~  0.066,  £2  ~  0.20,  and  f  ~  0.15.  The  larger  rotational  constants  of  CH3 
compared  to  C2H5  are  the  reason  for  larger  £2  and  £  in  CgHi0+  compared  to 
C9H12+.  In  spite  of  the  fact  that  these  estimates  are  only  approximate,  it  is  clear 
that  §i,  £2,  and  §  are  below  unity.  Therefore,  for  the  molecular  ions  considered 
here,  there  will  be  marked  deviations  in  the  dynamics  from  the  adiabatic  limit 
(reductions  of  kcap  up  to  about  40%,  see  [1],  [4],  and  [5]). 

Details  of  our  CT  calculations  have  been  described  previously  [1-6]  and 
are  not  repeated  here.  The  initial  conditions  of  the  trajectories  were  selected 
randomly  from  a  uniform  phase  space  distribution  compatible  with  the  respec¬ 
tive  E  and  J.  About  104  trajectories  were  required  to  reduce  statistical  scatter 
of  the  results  for  a  given  (E,  J)  to  less  than  2%.  Thermal  capture  rates  were 
determined  either  by  thermal  averaging  of  the  (E,  J)-results  or  by  thermally 
weighted  selection  of  initial  conditions  [1],  Both  methods  gave  identical  re¬ 
sults,  confirming  the  validity  of  the  numerical  procedure. 

Our  determination  of  the  number  of  open  channels  W(E,  J)  for  the  tran¬ 
sitional  modes  starts  from  the  well  known  state-counting  for  free  rotor-pairs 
of  joint  total  E  and  J28-30.  In  the  absence  of  centrifugal  barriers  this  defines 
W0 (E,  J).  The  result  for  large  J  is 


W0  (E,  large  J)  =  nE2/6ata2  (3.2) 

/=! 

where  Bu  and  B2i  are  the  rotational  constants  (in  energy  units)  and  cr,  and  cr2 
the  symmetry  numbers  of  fragments  1  and  2,  respectively.  For  small  J,  treating 
the  fragments  as  spherical  tops  (with  rotational  constants  Bx  and  B2),  one  has 

W0(E,  sma\\J)^8(2J+l)E5,2/l5ala2B,B2(B]+B2y/2 .  (3.3) 

State-counting  for  asymmetric  top  fragments  is  also  easily  done.  In  prac¬ 
tice,  an  empirical  interpolation  between  Eqs.  (3.2)  and  (3.3),  such  as  de¬ 
scribed  in  [28],  often  suffices:  with  y  —  W0(E,  J)/W0(E,  large  J)  and  x  = 
W0(E,  small  J)/W0(E,  large  J),  one  may  use 

y  ~  {0.08x1/2  exp[— l.l(x  -  l)2]  + 1}  tanhx .  (3.4) 

An  alternative  interpolation  between  Eqs.  (3.2)  and  (3.3)  using  erf(x)  was 
proposed  in  [31]. 

The  appearance  of  centrifugal  barriers  E0(J)  on  the  MEP  potential,  such  as 
given  by  Eqs.  (2. 10) — (2.12),  leads  to  the  PST  expression 

Wfsr(E,  J )  =  wpst(E,  J)  W0(E ,  J) . 


(3.5) 
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Fig.  3.  Capture  probabilities  for  an  isotropic  potential  (PST)  and  an  anisotropic  potential 
wrsT(E,  J )  (dashed)  and  w(E,  J )  (full),  respectively  (from  top  to  bottom  curves  for,  J  = 
25,  55,  75,  105,  135,  and  180  with  E0(.J)/hc  =  0.83,  12.8,  35,  99.5,  211,  and  481  cm-1, 
respectively;  results  for  C8H|0+  with  dipole-dipole  type  anisotropy  with  C2/D  «s  20). 


The  capture  probability  wfsT(E,  J)  was  either  determined  by  CT  calcula¬ 
tions  on  the  radial  potential  of  Sect.  2  or  by  state-counting  with  the  condition 
E  >  E0(J).  Figs.  3  and  4  include  the  results  for  the  CgH10+  system  for  the 
two  choices  of  anisotropy  used.  Like  in  our  earlier  work  [6],  for  two  molecular 
fragments, 

wkt(E,  J)~[l-E0(J)/E]2  (3.6) 

provides  an  adequate  approximation  (for  an  atom  and  one  molecular  fragment, 
the  exponent  2  in  Eq.  (3.6)  would  be  replaced  by  unity;  one  should  note  that  the 
centrifugal  barriers  E 0  depend  on  J  as  well  as  on  the  orbital  angular  momen¬ 
tum  with  quantum  number  L;  the  results  shown  for  w?ST(E,  J )  in  Figs.  3  and  4 
are  for  randomly  selected  L  and,  therefore,  are  slightly  scattered  and  somewhat 
below  [1  —  E0(J)/E]2). 

While  Eqs.  (3.2)-(3.6)  provide  easy  access  to  the  PST  results  for  WTST 
( E ,  J),  the  additional  constraints  introduced  by  the  anisotropy  of  the  po¬ 
tential  cannot  be  represented  in  such  a  simple  manner.  Fig.  3  compares  the 
capture  probabilities  in  PST,  w?ST,  with  capture  probabilities  w(E,  J )  for  the 
anisotropic  potential  such  as  derived  by  the  CT  calculations.  The  C8F11o+  sys- 
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tern  serves  as  the  example.  One  observes  that,  for  J  ->•  0  and  E  ->  E0(J), 
w(E,  J )  approaches  w?sr(E,  J).  In  contrast,  for  increasing  E  and  ./,  it  falls 
increasingly  below  wPST(£,  J).  In  the  language  of  transition  state  theory,  this 
would  correspond  to  transition  state  shifting:  with  increasing  excitation  en¬ 
ergy,  the  effective  bottle-necks  of  the  dynamics  smoothly  move  to  smaller 
interfragment  distances  where  the  potential  is  more  anisotropic  [32].  Conse¬ 
quently,  w(E,  J )  falls  below  w?ST(E,  J).  An  analogous  situation  applies  to 
angular  momentum:  with  increasing  J,  the  effective  bottle-necks  of  the  dy¬ 
namics  smoothly  move  into  a  more  anisotropic  range  of  the  potential.  Whereas 
Fig.  3  illustrates  w(E,  J )  for  dipole-dipole  type  anisotropy,  there  is  practically 
no  change  if  ion-dipole  type  anisotropy  is  chosen  instead,  as  shown  by  compar¬ 
ing  Figs.  3  and  4  (however,  the  parameters  C,  and  C2  differ).  The  results  for  the 
C9Hi2+  system  look  quite  similar  and,  therefore,  are  not  shown  here. 

An  alternative  representation  of  the  consequences  of  the  anisotropy  of  the 
PES  involves  “specific  rigidity  factors”  frisii(E,  J )  defined  by 

/rigid(£,  j)  =  vu(E,  J)/wrsT(E,  J ) .  (3.7) 

Fig.  5  demonstrates  fn$&{E,  J)  for  the  C8H10+  system.  With  increasing 
energy,  /rigid (E,  J )  falls  increasingly  below  the  PST  limit  (/pgsid  =  1).  Since 
W(E,  J)  for  the  transitional  modes  is  later  convoluted  with  contributions  for 
conserved  modes,  only  small  excitations  per  mode  are  of  importance  for  the 
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Fig.  5.  Specific  rigidity  factors  /^(is,  J)  for  the  CsH|0+-system  (curves  from  left  to 
right  are  for  J  —  25,  55,  75,  105,  135,  and  180,  respectively;  symbols:  results  from  CT 
calculations  for  dipole-dipole  type  anisotropy  with  C2/D  «  20;  full  curves:  analytical  ap¬ 
proximation  by  Eq.  (3.8)  with  m  =  2,  n  =  —2/3,  and  Essw/hc  140  cm~'). 


polyatomic  ions  considered  here.  Therefore,  the  low  energy  limiting  range  of 
/rigid  (E,  J )  is  most  relevant. 

For  convenience,  we  have  analytically  represented  the  derived  /rigid(Zi,  J ) 
in  the  form 


ME,  J)  ~  [1  +  —  E0(J)]/Esw}mY  .  (3.8) 

It  appears  worth  mentioning  that  the  calculated  rigidity  factors,  except 
for  small  60),  only  depend  on  E  —  E0(J)  and  do  not  show  additional  J- 
dependences.  In  addition,  the  exponents  m  and  n  were  found  to  be  close  to 
m  ~  2  and  n  ~  —2/3.  Keeping  these  values  of  m  and  n  fixed,  the  only  fit  pa¬ 
rameter  is  £sw.  Performing  the  CT  calculations  for  C8H10+  with  dipole-dipole 
anisotropy,  Esw/hc  ~  140  cm-1  was  fitted  which  corresponds  to  C2/D  —  20. 
Fig.  5  illustrates  the  agreement  of  the  CT  results  and  the  corresponding  empir¬ 
ical  fits  from  Eq.  (3.8).  Fig.  5  represents  the  connection  between  the  potential, 
e.g.  characterized  by  dipole-dipole  anisotropy  and  its  amplitude  factor  C2/D, 
and  the  rigidity  factor  characterized  by  E sw.  We  found  that  E sw  is  inversely 
proportional  to  C2/D.  Empirical  fits  of  Eq.  (3.8)  and  its  parameter  £sw  to  ex¬ 
perimental  results  then  also  lead  back  to  the  underlying  potential  parameter 
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C2/D  via  C2/D  ~  20  x  140  hccm~' / E^,.  Rigidity  factors  for  C9H12+  were 
found  to  look  quite  similar  to  those  shown  in  Fig.  5.  The  corresponding  relation 
for  C2/D  was  found  to  be  C2/D  «  18  x  1 10 /tc  cm_1/^sw- 

For  the  large  excitations  considered  in  the  present  work,  the  classical 
treatment  appears  adequate.  For  smaller  excitations,  individual  adiabatic  and 
axially-nonadiabatic  potential  curves  (including  Coriolis  coupling)  have  to  be 
determined  as  illustrated  in  [6]  and  [10].  This  difficult  task  is  beyond  the 
scope  of  the  present  article.  Instead  we  remain  within  the  SACM/CT  con¬ 
cept  and  convolute  W(E,  J)  from  the  CT  calculations  with  the  contributions 
from  quantized  conserved  oscillators.  It  appears  most  practical  to  do  this  by 
the  use  of  the  Beyer-Swinehart  algorithm,  using  the  derived  analytical  ex¬ 
pressions  for  W(E,  J)  of  the  transitional  modes  for  the  construction  of  the 
starting  array  in  the  algorithm  such  as  explained  in  Appendix  A  of  [33].  The 
frequencies  used  for  the  conserved  modes  of  C8H10+  (or  C9H12+),  corres¬ 
ponding  to  the  vibrations  of  C7H7+  and  CH3  (or  C2H5),  are  listed  in  Ap¬ 
pendix  A. 


4.  Specific  rate  constants  k(E,  J) 

Having  calculated  W(E,  J)  and  relying  on  the  usual  statistical  assumption  [8] 
of  one-channel  rate  constants  being  equal  to  1  /hp(E,  J),  the  SACM/CT  treat¬ 
ment  leads  to  specific  rate  constants  k(E ,  J)  =  Wl0t(E,  J)  x  [1  /hp(E,  J)]  as 
illustrated  below.  Rovibrational  densities  of  states  p(E,  J)  were  calculated  in 
the  harmonic  oscillator-rigid  rotor  approximation  using  fundamental  frequen¬ 
cies,  since  it  presently  is  not  feasible  to  properly  account  for  anharmonicity  in 
systems  of  this  size.  (Even  for  triatomic  systems,  apart  from  few  exceptions, 
see  [34],  anharmonic  rovibratonal  densities  of  states  near  to  the  dissociation 
limit  have  not  yet  been  determined  rigorously).  Vibrational  densities  of  states 
p{E,  J  —  0)  were  determined  by  Beyer-Swinehart  counting;  rotational  factors 
were  approximated  by  the  analytical  expressions  elaborated  in  [30]. 

The  resulting  k(E,  J)  for  the  dissociation  CSH10+  — >  C7H7+  +  CH3  are 
shown  in  Fig.  6.  The  figure  compares  PST  calculations  of  k(E,  J )  with  cal¬ 
culations  for  the  anisotropic  potential.  One  realizes  a  marked  influence  of  the 
“rigidity  of  the  potential”  due  to  its  anisotropy.  The  reduction  of  k{E,  J)  to 
values  below  k?SJ(E,  J )  becomes  increasingly  pronounced  with  increasing  en¬ 
ergy:  the  discrepancy  is  already  a  factor  of  20-30  at  E/ he  —  15  000  cm-1  and 
increases  to  almost  three  orders  of  magnitude  above  E/hc  =  40000  cm-1  (one 
should  note  that,  for  Figs.  6,  8,  and  10,  the  zero  of  the  energy  scale  is  the 
rovibrational  ground  state  of  the  molecular  ion  unlike  the  other  figures  of 
this  article).  The  increasing  drop  of  k(E,  J)  below  krsT(E,  J )  with  increas¬ 
ing  energy  is  the  sign  of  a  shifting  of  the  effective  dynamical  bottle-neck  of 
the  reaction  towards  smaller  interfragment  distances  where  the  anisotropy  is 
increasingly  more  pronounced.  For  a  fixed  energy  E,  there  is  only  a  weak  de- 
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Fig.  6.  Specific  rate  constants  k(E,  J )  for  the  dissociation  CsH,o+  ->  C7H7+  +  CH3  (up¬ 
per  full  curves:  k?ST(E,J)  for  J  =  0,  65,  and  150  from  bottom  to  top;  lower  dashed 
curves:  k(E,  J)  for  J  —  0  (middle),  65  (bottom),  and  150  (top),  see  text;  for  a  more 
detailed  representation  of  the  ./-dependence,  see  Fig.  7;  full  circle:  experimental  re¬ 
sult  from  [37];  full  curve  between  E  =  16000  and  22000cm"1:  experimental  results 
from  [38];  full  curve  from  k(E ,  / )  =  102  to  10s  s"1:  empirical  fitting  of  experimental  re¬ 
sults  by  Eq.  (6.2),  see  text). 


pendence  of  k(E,  /)  on  /.  This  is  shown  in  Fig.  6  and  illustrated  on  a  larger 
scale  in  Fig.  7.  The  weakness  of  the  ./-dependence  corresponds  to  that  cal¬ 
culated  for  the  fragmentations  of  neutral  toluene  and  ethylbenzene  [35]  and 
also  to  the  experimental  results  obtained  for  the  fragmentation  of  benzene 
cations  [36].  kPST(E,  J )  shows  the  same  ./-dependence  as  k(E,  J )  such  that  the 
./-dependence  of  fn&&(E,  J )  is  not  responsible  for  the  calculated  trends.  Since 
three  unrelated  factors,  i.e.  the  /-dependences  of  W¥sr(E,  /),  of  p(E,  J),  and 
of  E0(J),  govern  the  /-dependence  of  k(E,  J),  there  is  no  simple  explanation 
for  the  resulting  /-dependence. 

The  model  potentials  discussed  in  Sect.  2  contain  a  number  of  quantities 
that  needed  to  be  fixed  before  the  CT  calculations  were  made.  Apart  from  the 
known  molecular  parameters  given  in  Appendix  A,  the  Morse  parameters  ft  as 
well  as  the  anisotropy  amplitudes  CxjD  or  C2/D  had  to  be  specified.  Whereas 
/3  could  be  related  to  a  vibration  defining  the  reaction  coordinate,  C,  / D  or 
C2/ D  were  parameters  fitted  through  a  comparison  of  the  calculated  k(E,  /) 
with  experimental  results.  As  a  first  guess  of  the  order  of  magnitude  we  have 
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Fig.  7.  ./-dependence  of  k(E,  J)  and  k?ST(E,  J)  for  the  dissociations  C8H,0+  -*■  C7H7+  + 
CH3  (full  lines)  and  C9H|2+  ->  C7H7+  +  C2H5  (dashed  lines)  (modelling  results  for  PST 
(upper  curves)  and  for  anisotropic  potential  (lower  curves)  at  E/hc  =  30000  cm"1). 


estimated  these  parameters  by  means  of  the  anisotropy  parameters  from  the 
Appendix  of  [6].  In  this  way,  we  choose  C2/D  «  20  for  C8H10+.  This  value 
was  used  for  the  illustrative  examples  of  the  CT  calculations  in  Figs.  3-5. 
However,  this  value  of  C2/D  gives  k(E,  J )  which  differ  from  the  measure¬ 
ments.  In  order  to  obtain  a  more  realistic  C2/D,  the  corresponding  ESVJ  in 
Eq.  (3.8)  was  varied  until  an  optimum  experimental  value  of  k(E )  was  repro¬ 
duced.  Choosing  the  experimental  k(E )  near  32000  cm"1  from  [37]  for  this 
purpose,  the  fit  parameters  Esw/hc  ~  38  cm-1  and  C2/D  ~  74  are  obtained. 
In  this  way  the  complete  k(E,  J)  is  predicted  as  shown  in  Figs.  6  and  7.  The 
low  energy  results  near  18  000  cm"1  from  [38]  are  not  well  reproduced,  predict¬ 
ing  a  slightly  different  energy  dependence  of  k(E).  However,  the  differences 
still  appear  to  be  within  the  uncertainties  of  both  the  experiments  and  the 
model  potential.  Better  agreement  would  be  obtained  by  slightly  lowering  the 
value  of  E0  which  is  not  well  known.  In  addition,  the  disagreement  is  in  the 
range  of  k{E)  where  radiative  cooling  may  be  important.  This  is  not  modelled 
here. 

Modelled  k(E,  J )  for  the  dissociation  C9HI2+  ->  C7H7+  +  C2H5  are  com¬ 
pared  with  measured  k(E)  from  [39]  in  Fig.  8.  In  this  case,  dipole-dipole 
anisotropies  from  Eq.  (2.16)  were  used  in  CT  calculations  with  starting  values 
of  C2/D—  18  which  were  found  to  correspond  to  ESw/hc  =  110  cm"1  in 
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Fig.  8.  As  Fig.  6,  for  the  dissociation  C9H,2+  ->•  C7H7+  +  C2H5  (upper  dashed  curves: 
k?ST(E,  J )  for  J  =  0,  100,  and  200  from  bottom  to  top;  lower  full  curves:  k(E,  J)  for 
/  =  0  (middle),  100  (bottom),  and  200  (top);  see  text;  full  circles:  experimental  results 
from  [39]). 


Eq.  (3.8).  The  fit  to  the  experiments  then  was  done  choosing  the  experi¬ 
mental  point  from  [39]  near  32000  cm-1  which  fixes  ESw/hc&  17.3  cm-1. 
This  value  again  is  smaller  than  the  starting  value  of  our  CT  calculations 
and  corresponds  to  C2/D^  114.  The  complete  k(E,J)  is  constructed,  in 
this  case  agreeing  very  well  with  the  low  energy  point  near  18000  cm-1. 
The  modelled  /-dependences  of  k(E,  J )  show  the  same  trends  as  observed 
for  C8H10+  as  shown  in  Fig.  7.  The  shape  of  the  /-dependence  of  k(E,  J) 
slightly  differs  from  that  obtained  for  C8H10+.  The  details  of  this  trend  are 
again  caused  by  a  complicated  interplay  of  p(E,  /),  of  WrsT(E,  /),  and 
of  E0(J).  More  practical  aspects  of  the  modelled  k(E,  /)  are  discussed 
below. 


5.  Thermal  capture  rate  constants  kcap(T) 

The  number  of  open  channels  W(E,  J)  also  forms  the  basis  for  the  calculation 
of  thermal  capture  and  high  pressure  association  rate  constants.  For  example, 
the  PST  expression  of  Eq.  (2.14)  with  the  centrifugal  pseudo-partition  function 
of  Eq.  (2. 1 3)  and  the  centrifugal  barriers  from  Eq.  (2.11)  leads  to  the  Langevin 
expression  krc^  =kL  —  2jtq{a/ji)l/2  for  ion-induced  dipole  capture,  see  above. 


Classical  Trajectory  and  Statistical  Adiabatic  Channel  Study  . . . 


733 


For  the  present  short-range  valence/long-range  induced  dipole  potentials, 
a  number  of  changes  arise:  (i)  the  PST  capture  rate  constant  from  Eq.  (2.14) 
has  to  be  calculated  with  E0(J)  from  Eq.  (2.10)  instead  of  Eq.  (2.11),  (ii)  the 
general  capture  rate  constants  kcap(T)  have  to  account  for  the  specific  rigidity 
factors  of  Eq.  (3.8),  (iii)  one  has  to  check  whether  the  frequencies  of  the  con¬ 
served  modes  can  be  identified  with  the  frequencies  of  the  separated  fragments, 
and  (iv)  one  has  to  relate  kca{l(T)  with  experimental  high  pressure  association 
rate  constants  krcccc(T)  by  inspecting  which  fraction  of  the  electronic  states  of 
the  separated  fragments  leads  to  the  bound  adduct  under  consideration.  In  the 
simplest  case,  the  latter  relation  is  established  by 


*rec,oo  (T)=kQap(T)QM)/[Q«(B)Q*(C)],  (5.1) 


where  QA  are  the  respective  electronic  partition  functions  of  A  &  B  +  C.  As¬ 
sumption  (iii)  often  is  true  to  a  sufficient  extent,  because  the  major  part  of 
changes  of  conserved  mode  frequencies  from  A  to  B  +  C  will  generally  occur 
at  shorter  distances  than  relevant  for  the  effective  bottle-necks  of  the  transi¬ 
tional  mode  dynamics,  see  the  original  version  of  S ACM  [13].  In  this  case,  the 
calculation  of  kcap  will  only  concern  transitional  modes,  i.e. 

.  ,„_w(  h 2  y/2  q* 

CapC  j  h\2nnkTj  (2 rot  (B)  £> rot  (C)  {  ) 

with  the  pseudo-partition  function  Q*  given  by 


Q* 


»  7 

£(2J+1)  / 

J= 0  „ 


E  o(J) 


W  (E,  J)  exp  i-E/kT)  dE/kT 


(5.3) 


where  W{E,  J)  here  is  the  number  of  open  channels  of  the  transitional  modes 
and  (2  rat  are  the  rotational  partition  functions  of  the  separated  fragments  B 
and  C,  respectively. 

Evaluating  kpaspT(T)  for  the  capture  C7H7++  CH3  ->  C8H10+ withEqs.  (2.10), 
(2.13),  and  (2.14)  leads  to  the  results  shown  in  Fig.  9.  k?c^(T)  increases 
slightly  from  the  Langevin  rate  constant  at  low  temperatures  to  larger  values 
at  higher  temperatures  which  reflects  the  stronger  attraction  of  the  short- 
range/long-range  switching  potential  of  Sect.  2. 1  in  comparison  to  a  pure 
ion-induced  dipole  potential.  Including  the  specific  rigidity  factors  fTigu(E,  J ) 
from  Eq.  (3.8)  in  Eqs.  (5.2)  and  (5.3)  leads  to  the  curve  for  kcap(r)  in  Fig.  9 
which  falls  increasingly  below  kpaspT  (T)  with  increasing  temperature.  Again  this 
reflects  a  shifting  of  the  effective  bottle-necks  of  the  transitional  mode  dynam¬ 
ics  into  the  more  anisotropic  short-range  valence  portion  of  the  PES.  It  should 
be  emphasized  that  Fig.  9  shows  thermal  capture  rate  constants  which  are  fully 
consistent  with  the  specific  rate  constants  k{E,  J)  from  Fig.  6. 
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Fig.  9.  Thermal  capture  rate  constants  kap(T)  for  the  reaction  C7H7+  +  CH3  ->  CSH,0+ 
(upper  curve:  PST,  CT  calculations  and  Eqs.  (2.10),  (2.13),  and  (2.14);  lower  curve:  kav, 
CT  calculations  and  calculations  with  Eq.  (3.8)  fitted  to  reproduce  k(E,  J)  from  Fig.  8,  see 
text;  empirical  representation  by  Eq.  (7.2)  with  T0  =  30.4  K;  •:  experimental  points  from 
dissociation  experiments  converted  by  the  equilibrium  constant  from  [43]). 


6.  Simplified  modelling  of  k(E,  J) 

The  above  treatment  of  the  dynamics  on  the  level  of  the  SACM/CT  con¬ 
cept  provides  an  accurate  connection  between  k(E,  J)  and  the  PES.  In  the 
following  we  also  compare  these  results  with  a  series  of  simpler  procedures 
which  may  be  of  practical  use.  Since  one  has  to  fit  potential  parameters  in  any 
case,  one  may  also  ignore  the  potential  from  the  beginning  and  fit  k(E,  J)  to 
a  limited  number  of  experimental  points.  This  appears  desirable  when  k(E ) 
has  to  be  extrapolated  and,  in  particular,  when  kinetic  shifts  have  to  be  deter¬ 
mined  in  order  to  derive  bond  energies  E0  from  k(E).  It  was  suggested  in  [40] 
that  k(E),  for  polyatomic  dissociating  species,  can  be  approximated  by  an 
expression 

k(E)*A*(E-E0f-'  (6.1) 

with  the  true  bond  energy  E0,  an  empirical  exponent  s  *  —  1  and  a  constant 
factor  A*.  The  form  of  Eq.  (6.1)  should  be  distinguished  from  classical  Kas¬ 
sel  theory  which  has  an  additional  energy  dependence  in  the  denominator  and 
identifies  s*  with  the  number  of  activated  complex  oscillators.  After  rewriting 
Eq.  (6.1)  in  the  form 

k(E)  ~  A[(E-E0)/(ETzf-E0)Y'-'  , 


(6.2) 
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we  find,  for  C8H10+  ->■  C7H7+  +  CH3,  that  k(E)  could  well  be  represented 
with  A  =  1.6  x  108s-1,  Ea/hc  —  14070  cm-1,  Eref/hc  =  32420  cm-1,  and 
s*  =  5.962,  using  the  measured  k(E)  at  Ek(  from  [37]  as  the  reference  fit¬ 
ting  point  and  using  the  experimental  results  from  [38]  for  fixing  s*.  As 
seen  in  Fig.  6,  a  slightly  steeper  energy  dependence  of  k(E)  over  the  range 
16000-21 000  cm-1  in  [38]  was  determined  experimentally  than  given  by 
Eq.  (6.2).  While  the  fully  empirical  Eqs.  (6.1)  and  (6.2)  do  a  fair  job  in  re¬ 
producing  measured  and  calculated  k(E,  J  —  0),  the  ./-dependence  of  k(E,  J) 
escapes  a  simple  empirical  representation.  In  any  case,  the  weak  ./-dependence 
illustrated  in  Fig.  7  indicates  that  the  measured  k(E)  from  [37]  and  [38]  are  not 
different  from  k(E,  J  =  0)  within  experimental  error. 

Comparison  of  the  simplified  SACM  calculations  of  k(E,  J  —  0)  for  the 
dissociations  of  neutral  toluene  to  benzyl  +  H  or  phenyl  +  CH3  and  of  ethyl¬ 
benzene  to  benzyl  +  CH3  [35]  to  results  obtained  with  pure  model  valence  po¬ 
tentials  confirms  the  good  performance  of  Eq.  (6.1).  In  those  cases,  ./-depend¬ 
ences  of  k(E,  J )  were  determined  with  increasing  k(E,  J )  for  increasing  J  in 
C-C  bond  breaking,  like  in  the  present  work,  and  decreasing  k(E,  J)  in  C-H 
bond  breaking  [35],  in  contrast  to  the  present  work. 

Reliable  experimental  values  of  k(E)  at  three  energies  are  enough  to  de¬ 
termine  the  three  empirical  parameters  in  Eq.  (6.1)  and  therefore  sufficient  to 
derive  E0.  For  a  less  empirical  approach  and  for  a  determination  of  E0,  two 
sufficiently  distant  experimental  points  of  k(E)  would  suffice  if  k(E,  J)  from 
a  PST  calculation  is  combined  with  /ngld(£,  J )  from  Eq.  (3.8)  and  the  two  pa¬ 
rameters  E0  and  Es w  are  fitted.  While  the  full  SACM/CT  treatment  requires 
large  computing  efforts,  this  simple  “modified  PST  treatment”  is  by  far  less 
time  consuming.  Including  centrifugal  barriers  E0(J)  as  defined  in  Eq.  (2.10) 
into  the  PST  treatment,  allows  the  correct  ./-dependence  of  k(E,  J )  to  also  be 
satisfactorily  accounted  for. 

Common  practice  is  to  model  k{E,  J  =  0)  with  RRKM  theory  using 
oscillator-type  activated  complexes,  see  e.g.  [8,37,39],  Because  of  the  large 
number  of  adjustable  frequencies,  this  approach  can  obviously  also  reproduce 
measured  k(E).  However,  there  is  no  rigorous  justification  for  this  procedure 
in  simple  bond  fissions  without  reverse  barriers,  and  the  chosen  frequencies  of 
the  transitional  modes  have  no  direct  relation  to  the  properties  of  the  potential. 
Fixing  the  frequencies  of  the  conserved  modes  to  those  of  the  fragment  modes, 
one  might  try  in  RRKM  theory  to  treat  the  transitional  modes  as  rotors,  such  as 
done  in  PST.  However  this  causes  apparent  increases  of  the  rotational  constants 
in  order  to  mimic  rigidity  [41],  Fitting  this  increase  would  lead  to  an  energy  in¬ 
dependent  rigidity  factor  which  obviously  does  not  correspond  to  the  results  for 
k{E,  J)  shown  in  Figs.  6  and  8.  A  microcanonical  VRC-VTST-treatment  [18] 
of  the  transitional  modes  (see  above),  of  course,  would  produce  better  results. 
However,  it  also  requires  a  full  reduced-dimensionality  PES,  is  fairly  involved, 
requires  dynamical  corrections  [19],  and  needs  parameter  fitting  as  long  as  no 
ab  initio  potential  is  available. 


736 


J.  Troe  et  al. 


Fig.  10.  Simplified  representations  of  specific  rate  constants  k(E,  J  =  0)  for  the  dissoci¬ 
ation  C9H,2+  ->  C7H7+  +  C2H5  (see  Fig.  8,  experimental  points  from  [39];  heavy  full 
line:  empirical  fit  by  Eq.  (6.2)  with  A  —  3.6  x  10s  s"1,  E0/hc  —  13950  cm"1,  ErtS/hc  = 
32260  cm"1  and  s*  =  7.541;  light  full  line:  RRKM  fit  from  [39];  dotted  line:  SACM/CT 
modelling  adjusted  to  the  experimental  point  near  32000  cm"1  and  represented  by  modi¬ 
fied  PST  calculation  with  Eq.  (3.8)  with  ESw/hc  =  17.3  cm"1;  dashed  line:  PST  calcula¬ 
tion). 


Fig.  10  compares  k(E,  J  =  0)  from  the  various  approaches  for  the  reac¬ 
tion  of  C9H12+  — ►  C7H7+  +  C2  Hs.  Included  in  this  figure  are  an  empirical 
fit  using  Eq.  (6.2),  a  modified  PST  modelling  using  Eq.  (3.8)  being  equivalent 
to  the  full  SACM/CT  results,  an  RRKM  fit  from  [39],  and  the  experimental 
points  from  [39].  The  differences  of  the  different  approaches  become  clearly 
visible  only  when  larger  energy  ranges  are  considered.  However,  larger  energy 
ranges  generally  are  relevant  when  kinetic  shifts  of  specific  rate  constants  are 
estimated  [8],  Therefore,  the  different  simplified  methods  described  in  this  sec¬ 
tion  may  also  lead  to  different  results  among  each  others  and  differ  from  the 
detailed  SACM/CT  results. 


7.  Simplified  modelling  of  kcap(T) 

After  having  presented  simplified  ways  to  characterize  k(E,  J),  one  may  also 
try  to  represent  kcaf(T)  in  a  simple  manner  that  is  consistent  with  the  simpli¬ 
fied  k(E,  J).  For  the  C8H10+  system,  the  decrease  with  increasing  tempera- 
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Fig.  11.  As  Fig.  9,  but  for  C7H7+  +  C2H5  ->  C9HI2+  (upper  curve:  PST;  lower  curve: 
ka p  corresponding  to  k(E,  J)  from  Fig.  8;  experimental  point  at  603  K  from  high  pres¬ 
sure  dissociation  experiments  of  C9H12+  after  conversion  by  the  equilibrium  constant  [42]; 
empirical  representation  by  Eq.  (7.2)  with  T0  =  13.8  K;  minor  changes  of  the  employed 
equilibrium  constant  would  bring  the  point  and  Eq.  (7.2)  into  agreement,  see  [42]). 

tures  of  the  capture  rate  constants  from  values  dominated  by  the  long-range 
electrostatic  potential  to  values  dominated  by  the  short-range  valence  com¬ 
ponent  of  the  potential  is  shown  in  Fig.  9.  Similar  results  were  obtained  for 
the  C9H12+  system  and  are  illustrated  in  Fig.  1 1.  As  a  result  of  the  anisotropy 
of  the  potential,  which  is  calibrated  by  the  measurements  of  k(E),  &cap  falls 
markedly  below  the  results  from  PST  in  both  cases.  Before  asking  whether 
this  prediction  is  confirmed  by  experiments,  we  try  to  provide  an  analytical 
expression  for  kQap(T)  which  is  linked  to  the  values  of  E sw  characterizing 
k(E ,  J).  Calculating  kcv(T)  for  the  short-range  potential  separately  would  lead 
to  k™(T).  One  may  expect  that  this  is  the  limit  of  kcap(T)  approached  at 
high  temperatures  while  k?0^(T)  is  the  limit  approached  at  low  T.  We  found 
that  the  transition  between  kpJJ(T)  and  k^/(T)  then  can  be  approximated 
by 


kcap  ( T )  *  k%  (T)  +  [^;  (T)  - k*  (T)]  [T0/  (To  +  T)\ L5  .  (7. 1) 

As  k^(T)  k™(T)  in  the  present  case,  Eq.  (7.1)  simplifies  to 

k^(T)^k™(T)  [T0/ (T0 +  T)]'  S  . 


(7.2) 
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In  contrast  to  the  value  1 .5  of  the  exponent,  we  found  empirically  that  the 
parameter  T0  is  linked  to  Esw  by 

kT0  ~  0.56£sw  (7.3) 

where  £sw  is  from  Eq.  (3.8).  E& w  is  determined  either  through  theoretical  mod¬ 
elling  or  through  empirical  fitting  of  k(E).  Eqs.  (7. 1)— (7.3)  then  provide  the 
internally  consistent  link  between  k{E)  and  kcap(T).  At  the  same  time,  through 
the  connection  between  T0,  £sw,  and  C2/D ,  the  link  to  the  dominant  parameter 
characterizing  the  anisotropy  of  the  potential  is  established. 

The  prediction  that  thermal  rate  constants  for  capture  of  CH3  or  C2H5  by 
C7H7+  are  far  below  the  Langevin  rate  constants  for  ion-induced  dipole  capture 
needs  experimental  verification.  The  corresponding  experiments  have  not  yet 
been  made  for  the  recombination  process.  However,  recent  turbulent  ion  flow 
tube  (TIFT)  measurements  of  the  rate  of  pyrolysis  of  C8H10+  and  of  C9H12+ 
have  been  reported  in  [42]  and  [43].  These  experiments  have  been  made  very 
near  to  the  high  pressure  limit  of  the  unimolecular  dissociation.  Converting 
these  results  by  the  equilibrium  constant  to  the  high  pressure  recombination 
rate  constant,  one  obtains  krec  oo  which  in  this  case  indeed  arc  close  to  kcap  as 
derived  from  the  measured  k(E) .  The  corresponding  results  are  included  in 
Figs.  9  and  1 1 .  The  near  agreement  confirms  the  internal  consistency  of  the 
present  approach,  i.e.  the  measured  thermal  dissociation  rates  are  consistent 
with  the  experimental  specific  rate  constants  and  both  quantities  can  be  related 
to  the  relevant  potential  parameters.  The  small  differences  may  be  explained  by 
experimental  uncertainties  either  of  k(E)  or  of  kiKi(T),  or  by  E0  which  enters 
the  equilibrium  constants. 


8.  Conclusions 

In  the  present  work  specific  rate  constants  k(E,  J )  and  thermal  capture  rate 
constants  kcap(T)  for  representative  examples  of  molecular  ion  dissociations 
and  the  reverse  associations  have  been  calculated  in  a  consistent  way  by 
the  SACM/CT  approach,  Model  potentials  of  mixed  short-range  valence  and 
long-range  electrostatic  character  were  employed  and  direct  links  between 
rate  constants  and  potential  parameters  were  established.  The  chosen  exam¬ 
ples  C8H10+  and  C9H12+  are  characterized  by  small  Massey  parameters  for 
the  transitional  modes  such  that  their  dynamics  is  nonadiabatic.  The  same  is 
true  for  a  large  series  of  other  neutral  or  ion  dissociations  and  recombina¬ 
tions  such  as  H  +  CH3  44CH4  (§  =  0.35),  H  +  C2H3  44C2H4  (£  =  0.12), 
H  +  C2H5  44  C2H6  (£  =  0.088)  or  2CH3  44  C2H6  (f  =  0.95),  2C2H5  44  C4H,0 
(i  =  0.36),  and  2C7H7  44  C14H14  ($  =  0.21),  NC  4-  NO  -►  NCNO  (£  =  0.43), 
C6H5+  +  H  -►  C6H6+  (£  =  0.053),  and  C6H5+  +  Br  -*  C6H5Br+  ($  =  0.20).  If 
these  reactions  were  treated  by  VTST,  i.e.  through  phase  space  sampling  along 
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c.o.m.  reaction  coordinates,  dynamical  corrections  for  nonadiabacity  would 
have  to  be  applied,  see  [19]  and  [44], 

As  ab  initio  calculations  of  the  potential  for  molecular  ions  of  the  size  con¬ 
sidered  here  are  only  rarely  available  with  adequate  precision,  some  potential 
parameters  in  general  have  to  be  fit.  In  this  situation  it  may  suffice  to  do  a  pa¬ 
rameter  fitting  directly  to  empirical  expressions  representing  the  experimental 
data.  We  have  proposed  several  simple  ways  which  minimize  the  number  of  fit 
parameters.  These  have  been  validated  by  comparison  with  the  full  dynamical 
SACM/CT  treatment.  k(E,  J)  and  kcap(T)  for  both  the  full  theoretical  treat¬ 
ment  and  the  simplified  parametrized  level  were  determined  in  a  consistent 
manner.  In  this  way,  the  consequences  of  having  a  short-range  valence  and 
long-range  electrostatic  character  of  the  potential  could  be  elaborated  in  detail. 
A  consequence  of  the  marked  anisotropy  of  the  short-range  valence  potential 
are  specific  rate  constants  k(E,  J)  which  are  far  below  the  predictions  from 
phase  space  theory.  Likewise,  thermal  recombinations  of  CH3  and  C2H5  with 
C7H7+  are  predicted  to  have  rate  constants  far  below  the  Langevin  rate  constant 
in  the  high  pressure  limit.  The  latter  prediction  was  confirmed  by  pyrolysis  ex¬ 
periments  of  C8H10+  and  C9H12+.  Finally,  the  SACM/CT  calculations  of  the 
present  work  have  been  complemented  by  product  energy  distributions  in  [22] 
which  are  internally  consistent  with  the  present  results  for  k(E,  J )  and  kcap(T). 
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Appendix:  Molecular  parameters 

a)  C8H10+  — >  C7H7+  +  CH3  (parameters  taken  from  [37]  and  [38]): 

E0/hc  =  14070  cm-1;  5,(C7H7+)//*c  =  0.18, 0.094,  and  0.062  cm-1;  ct(C7H7+) 
=  2;  f?;(CH3)//*c  =  9.3,  9.3,  and  4.65  cm-1;  cr(CH3)  =  6;  Morse  parameter 
P  =  hccoJn/2h2E0  =  3.68  A-1  with  a>  =  1003  cm-1;  re  -  3.46  A;  a(CH3)  = 
2.334  A3  from  [45];  frequencies  w,/cm3  =  163,  337,  342,  407,  510,  582,  611, 
628,  778,  782,  843,  946,  963,  976,  1000,  1012,  1046,  1089,  1112,  1134,  1157, 
1307,  1311,  1326,  1420,  1452,  1540,  1547,  1583,  2975,  3011,  3015,  3017, 
3036,  3039,  3067  for  C7H7+,  and  304,  1365,  1365,  2911,  3072,  3072  for  CH3; 
P re  =  12.73,  prBe  =  7.123,  /3rsw  =  14.186,  aq2p4/2D  =  1743. 
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b)  C9H12+->  C7H7+  +  C2H5  (parameters  taken  from  a)  and  from  [39]  and  [42]): 

E0/hc  —  13950  cm-1;  J?,(C2H5)/Ac  =  3.51,  0.75,  and  0.70;  <t(C2H5)  =  1; 
Morse  parameter  ^  =  5.13  A-1  with  co  =  1061  cm~!;  re  =  4.74  A;  a(C2H5)  = 
4.27  A3  from  [45];  frequencies  &>,/ cm3  =  155,  316, 784,  963,  997, 1166,  1383, 
1437,  1448,  1459,  2839,  2876,  2925,  2970,  3064  for  C2H5;  £r„  =  24.32, 
=  12.43,  /3rsw  =  26.34,  ctq2pA/2D  =  12  300. 
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